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Abstract. We study a class of holomorphic matrix models. The 
integrals are taken over middle dimensional cycles in the space 
of complex square matrices. As the size of the matrices tends to 
infinity, the distribution of eigenvalues is given by a measure with 
support on a collection of arcs in the complex planes. We show 
that the arcs are level sets of the imaginary part of a hyperelliptic 
integral connecting branch points. 



1. Introduction 

According to Dijkgraaf and Vafa [HE], the effective glueball superpo- 
tential of M = 1 supersymmetric U(n) gauge theory has an asymptotic 
expansion given by the planar part of the topological expansion of a 
matrix model. To give the effective potential for various vacua of the 
gauge theory, Dijkgraaf and Vafa propose a formula obtained by saddle 
point expansion of the matrix integral around different critical points. 
The filling fractions, i.e., the fraction of the eigenvalues sitting close to 
each of the critical points of the potential are the parameters selecting 
the vacua. 

In this paper we give a way to define non-perturbatively (i.e., beyond 
the saddle point expansion) the matrix integrals. Different integrals are 
integrals over different cycles in the space of normal matrices. In the 
case of generic polynomial potentials it is possible to construct a cycle 
for each of the critical points, so that all effective potentials consid- 
ered by Dijkgraaf and Vafa arise in the asymptotic expansion of our 
integrals. We note that the idea of integrating over eigenvalues along 
suitable contours in the complex plane appears in special situations in 
|3| . David's contours were also more recently considered in j3] , where a 
(non-holomorphic) modification of the integral was proposed to obtain 
arbitrary filling factors. 



This paper is based on the ETH diploma thesis of the second author. 
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A second result of this paper is the description of the asymptotic 
distribution of eigenvalues for each critical point. The eigenvalues, as 
predicted by Dijkgraaf and Vafa, lie asymptotically along arcs connect- 
ing branch points of a two-fold cover of the complex plane. We give a 
reality condition which specifies the shape of the arcs and the density 
of eigenvalues. 

Before considering the case of general NxN matrices, it is instructive 
to consider the case N — 1 of ordinary integrals. Consider the Airy 
integral 



This integral gives a formal solution of the Airy differential equation 
—p 2 u"(x) + xu(x) =0. Indeed the saddle point expansion around the 
two critical points of the integrand give the two linear independent 
formal power series solutions. The non-perturbative solutions can be 
obtained by integrating along paths in the complex plane connecting 
any two of the three directions at infinity where the integrand goes to 
zero. More generally, the integrals of the form J exp(^(xt + P(t))dt 
for a generic polynomial P(t) of degree n are solutions of a differential 
equation of order n — 1. We show that a system of n — 1 linearly 
independent solutions are obtained by choosing n integration contours 
in the complex plane going to infinity in directions where the integrand 
decays exponentially. Moreover the contours can be chosen so that each 
of them passes through exactly one of the n critical points. It follows 
that the n — 1 saddle point expansions at each of the critical points 
are the asymptotic expansions of true solutions given by convergent 
integrals. 

In the case of matrices we consider integrals of the form 



over real iV 2 -dimensional cycles V in the space of complex NxN 
matrices. The potential p is a polynomial with complex coefficients, 
g?$ = Aj t kd§jt and the normalization factor is such that (1)jv = 1- 
The observables F(<&) are holomorphic functions invariant under con- 
jugation. In the well-studied case of integrals over hermitian matrices, 
with a polynomial p with real coefficients bounded below, the saddle 
point approximation becomes exact in the large N limit with fixed p 
and the relevant critical point is described by an asymptotic density of 
eigenvalues which has support on a union of intervals on the real axis. 
If p is small, these intervals are small neighborhoods of the minima of 
the potential p. To obtain the densities of eigenvalues needed to make 
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contact to gauge theory one considers the variational problem for criti- 
cal points subject to the side condition that the fraction of eigenvalues 
in the vicinity of each of the critical points (not just minima) are given 
numbers. 

We consider here the case of a generic polynomial p of degree n 
with complex coefficients, with distinct critical points zi, . . . , z n -i, and 
propose to consider integrals over cycles in the space of normal matrices 
(a matrix is normal if it commutes with its adjoint or, equivalently, if 
it is conjugated to a diagonal matrix by a unitary matrix). The cycles 
we consider are parametrized by integers N 1: . . . , N n _i whose sum is N 
and are characterized by the condition that N k eigenvalues belong to a 
path Tfc in the complex plane going through z k and going to infinity in a 
direction where Re(p(z)) — > oo. In the limit N — > oo with = N^/N 
fixed, the eigenvalues in the saddle point approximation (supposed to 
be exact in the limit) are distributed along n — 1 arcs in the complex 
plane. For small fi the arcs are close to the critical points and the kth 
arcs contains the fraction of the eigenvalues. 

2. The one-dimensional case 

Let p(z) = a n z n + • • • + a±z + be a polynomial of degree n > 2 
with complex coefficients. We want to consider integrals of the form 

^q(z)e~l P ^dz, fi>0. 

for polynomials q(z). Before considering the question of integration 
cycles we may evaluate such integrals as asymptotic series as ji — > by 
formal application of the saddle point method at each of the critical 
points Zi, . . . , z n -i, which we assume to be distinct. In this way we get 
n—\ asymptotic series of the form exp(— p(zk) / /-f)(ci + c^^i + • • • ) and 
the question is whether these are asymptotic expansions of our integral 
for suitable cycles rV 

The cycles which we should consider here are (linear combinations 
of) paths for which the integral converges. As the integrand is holo- 
morphic, homotopic paths will give the same answer and what matters 
is the behavior at infinity. As z — > oo, p(z) ~ a n z n , so there are n 
directions in the complex plane for which Re(p(z)) — > +oo as z tends 
to infinity in these directions. Let us call these asymptotic directions 
valleys as in these directions the integrand decays exponentially. Neigh- 
boring valleys are separated by hills, which are directions of exponential 
increase of the integrand. So the cycles one needs to consider are linear 
combinations of infinite paths connecting pairs of distinct valleys. As 
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paths connecting two valleys can be deformed into sums of paths con- 
necting the two valleys with any third one, there are only n — 1 linearly 
independent cycles. 

Let us assume for simplicity that the critical values p(zk) have dis- 
tinct imaginary parts. Then there is a canonical way to associate to 
each critical point Zk a path Tk in such a way that the asymptotic ex- 
pansion of the integral over Tk as \i — > is obtained by the saddle point 
expansion at z^. Namely, we take the steepest descent paths(see, e.g., 
j3] and [U]) emerging from z^, defined by the condition that the tangent 
vector at each point points in the direction of the gradient of R,e(p(z)). 
As p(z) is holomorphic, the gradient of Rep(z) is orthogonal to the 
gradient of lmp(z) by the Cauchy-Riemann equations. Thus steepest 
descent paths are level lines for the imaginary part of p(z). Each non- 
degenerate critical point Zk is at the intersection of two such level lines. 
One of these two lines, the one along which ~Re(p(z)) takes its minimum 
at Zk, is the steepest descent path T k . Along the other line T' k , the real 
part takes its maximum at z^. We claim that Tk is a smooth path going 
to infinity in both directions and connecting two valleys. Indeed, Tk 
is (in suitable parametrization) given by a solution of the differential 
equation z(t) =p'(z(t)). It follows that f t Re(p(z(t))) = \p'(z(t))\ 2 . By 
our assumption, critical values have distinct imaginary parts, so the 
steepest descent path passing through z k may not come close to any 
other critical point. Thus \p'(z(t))\ is bounded below so that, as we go 
away from Zk, Rep(z(t) must go to infinity and the path Tk connects 
two valleys. Similarly T' k connects two hills. As Tk and T' k cross at 
Zk, the two valleys connected by Tk are separated by hills and are thus 
different. 

We have thus shown that for each critical point Zk there is a steepest 
descent path Tk going through Zk and connecting pairs of different 
valleys. On Tk the real part of p(z) is minimal at Zk so that the saddle 
point expansion at Zk indeed gives the asymptotic expansion of the 
integral over IV 



3. Matrix integrals 

3.1. Integration cycles. We consider matrix integrals of the form (Q) 
for p(z) a polynomial of degree n. They are integrals of holomorphic 
differential forms over N 2 dimensional cycles. For each set of natural 
numbers Ni, . . . , N n _i summing up to N, we have a cycle in the normal 
matrices, characterized by the condition that Nk eigenvalues run over 
the path T k of the previous section. More precisely, the cycle T is 
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parametrized by R N <g> U(N)/U(1) N : 

( Ai(ti) 

(t,U)»U ••. | 

The first iVi diagonal elements Ai(i), . . . , AjVi(^) are parametrizations 
of r 1; the next N 2 are parametrizations of T 2 and so on. 

The usual argument to reduce the integral to an integral over the 
eigenvalues (see 9j) gives 

(F)n = ^T f Ffri, • • • , Aiv)e~ f 52Up( X ^ J](A, - A fc ) 2 dAx • ■ • d\ N . 
Z x J j<k 



The integral is over T 1 1 x ■ ■ • x T n " 1 ^ 1 and the function F, a function 
on matrices invariant under conjugation, is regarded as a symmetric 
function of the eigenvalues. 

3.2. The loop equation. To study the large N limit, it is useful to 
introduce the trace of the resolvent: 

as products of such traces are generating functions of polynomial func- 
tions invariant under conjugation. The "loop equation" |H] for this 
quantity is 

fx(uj(z) 2 ) N -p'(z){u(z)) N + f N {z) = 0, 

where 

f (z] _ 1 / V P / (^)-p / (A i 

\ i IN 

is a polynomial of degree n — 2 with leading coefficient na n . This 
equation can be derived from the identity 




3>fc 



In the limit N — > oo, the matrix integral is supposed to be dominated 
by an integral over a region where the eigenvalues are close (for small 
H) to critical points of p. With our choice of integration cycles and 
keeping rij = Nj/N fixed as iV — > oo, there will be a fraction rij of the 
eigenvalues close to Zj. 
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In the limit iV — > oo one expects that the saddle point approximation 
becomes exact and thus (u(z))n converges to 

(u(z)) = I —^—dv(z'), 
J z z 

for some probability measure v with support in regions around the 
critical points of p and so that nk is the measure of the region around 
Zfr. Technically one assumes that the limit iV — > oo exists and that 
lim N ^ 00 ((uj(z)' 2 ) at — (u(z))%) = 0. The function (uj(z)) is defined and 
holomorphic for z outside the support of the measure. Setting 

(2) y(z)=2p(u(z))-p'(z), 
we finally obtain 

(3) y(z) 2 =p'(z) 2 -Apf(z), 

for some polynomial f(z) of degree n — 2 with leading coefficient na n . 
Thus the function y(z), which is a priori defined on the complement of 
the support of the measure, has an analytic continuation to a two-fold 
covering of the complex plane. The original function y is the branch of 
the function defined by (jSJ) which behaves at infinity as —p'(z). 



3.3. The density of eigenvalues. With the analogy with the case of 
hermitian matrices in mind, it is reasonable to assume that the measure 
v has support on a collection of arcs 71, ... , 7„_i which for small p are 
close to the critical points of p, and that du(t) = p(t)dt for some density 
of eigenvalues p defined on the arcs: namely, the measure u(U) of a set 
U intersecting one of the arcs, say 7,-, in a piece jj n U is 

u(U) = [ p(z)dz. 

Note that for the right-hand side to be defined we need to fix an ori- 
entation on 7j. Then y(z) is a holomorphic function outside on the 
complement of the arcs and the density on jj is related to the discon- 
tinuity of y: 

pO) = — j^- (y(x~) - y(x + )), x e t,-. 

Here y(x + ) (y(x~)) denotes the limit of y(z) as z tends to x from the 
left (from the right) of the oriented curve jj. 

From this information we deduce that the arcs jj connect pairs of 
zeros of y, the branch points of the hyperelliptic curve y 2 = p'(z) 2 — 
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4pf(z) of genus n — 2. The measure of the jth arc is then the period 
1 



(4) rij = I y(z)dz, j = 1, . . . , n - 1, 



4* 

over a cycle A,- enclosing the pair of branch points in counterclockwise 
direction. It follows from the condition on the leading coefficient of / 
that Yl n j = 1- We also note that since y(x~) = —y(x + ), we have the 
formula p{x) = y{x) jimp. 

3.4. The shape of the arcs. There remains to determine the precise 
form of the arcs and the coefficients in / as functions of the filling 
fractions rij. First of all, the relation between the n — 2 free complex 
coefficients 6 ,...,6„_3 of f(z) (recall that the leading coefficient is 
fixed to be 6 n _ 2 = na n ) and the periods rij (1 < j < n — 1) subject 
to ^2 rij = 1 is, locally around any point where the branch points are 
distinct, a holomorphic diffeomorphism, since the Jacobian matrix 



drij I f z 



v fc-i 



-dz, j, k — 1, . . . , n — 2. 



<% fc _i Airifx J A . y 

is non-degenerate, being the matrix of a-periods of a basis of Abelian 
differentials on a smooth curve. It follows that locally there exists a real 
(n— l)-dimensional submanifold in the complex space C n of coefficients 
bj which maps to real positive rij. 

The condition that fixes the shape of the arcs is the reality and 
positivity condition for the density: if 1 1— > z(t) is a parametrization of 
the arc 7j respecting its orientation, the condition is 

(5) P(*(*))^ > 0. 

Using = y(z~)/2nifi, we may parametrize the arcs (away from the 
endpoints) to be solutions of the differential equation 



z(t)=iy(z(t)-), 

connecting branch points. Alternatively, arcs may be described as level 
lines of a function: introduce the hyperelliptic integral 



i r 

F(z) = / y(z)dz. 

1 ; Amp J ZQ yK ' 



It is a holomorphic many-valued function on the complement of the 
support of the measure. As we go around an arc 7^, F(z) increases 
by rij so ImF(z) is single valued. The measure of a piece between 
two points xi, x 2 on a curve jj is F{x^) — F(x±) which is real. Thus 
the arcs jj are level lines of the imaginary part of F(z). Around a 
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branch point z which is a simple zero of y 2 , F(z) ~ const + (z — z ) 3 ^ 2 . 
Therefore there are three smooth level lines of F emerging from every 
simple branch points. In the most general situation the support of the 
measure may then be a graph consisting of level lines of ImF joining 
branch points. A more precise description is possible in the case of 
small fi to which we turn. 

3.5. Small coupling. For small /J > we claim that the arcs and the 
density of eigenvalues are determined completely by the filling fractions 
Tij through Q and (0). To show this, notice first that as fi — > 0, pairs 
of branch points z'j, z" converge to the critical points Zj and the periods 
7i j (eq. regarded as functions of \i and the coefficients bj of / are 
holomorphic at fi = 0. We have 

Since there is a bijective holomorphic correspondence between val- 
ues f{zj) at the distinct points zj and coefficients bj of /, we have 
at fi = 0, and by analyticity also for small /i, a biholomorphic map 
(b , . . . , b n _ 3 ) — > (ni, . . . , n n -i), n j — 1- m particular, we can invert 
this map and find a unique / for each set of rij > 0, such that Y2 n j = 1- 
It remains to show that for all small /i > there is a level line 7^ of 
ImF(z) connecting z'- to z". As fi — > 0, fiF(z) — ► — ^p(z) + const. 
The level lines of /xlmF(z) at /x = are thus the level lines of Kep(z). 
In the neighborhood of a non-degenerate critical point Zj they look like 
the left picture in Fig. ^ on any small circle around Zj each value is 
taken on at most four times. For positive small \i the critical point 





Figure 1 . Level lines of Im F for /i = (left) and for 
jj, > (right) 

Zj splits into two branch points z'j, z"- from each of which three level 

" J J 

lines emerge. The condition that the period rij is real implies that 
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ImF(z'j) = lmF(z"). The function g(z) = Im (F(z) — F(z'j)) is de- 
fined up to a sign in a neighborhood of z'j,z", so that its zero level 
line is uniquely defined. It follows that there is a level line jj of Im F, 
namely the zero set of g, joining z'j to z'j as on the right picture in 
Fig. ^ if none of the level lines emerging from z'j and z" were to join, 
g(z) would take the value zero six times on any small circle encircling 
z'p z'j, which cannot be, as this does not happen at [A — 0. Also, a level 
line cannot go from a point z'j or z'j to itself as the real part of F is 
monotonic along level lines. 

We conclude that for any small /i > 0, and any given filling fractions 
ni, . . . , n n _i > summing up to 1, there is a unique polynomial f(z) = 
na n z n ~ 2 + b n ^ 3 z n ~ 3 + ■ ■ ■ + bo, so that the curve y 2 = p'(z) — 4^f(z) has 
a-periods rij. The zeros of y are connected in pairs by arcs jj obeying 
the reality condition (j3j). According to the discussion above, these arcs 
are are the support of the measure and the density of eigenvalue is 
(2ni/j,)- 1 y(x-), x G jj. 

4. Concluding remarks 

We have given a non-perturbative definition of the matrix integrals 
that in the large N limit give the superpotentials considered in [HEj- 
We considered the case of a generic polynomial potential with complex 
coefficients. For small 't Hooft coupling /i, the density of eigenvalues 
was shown to be given by arcs connecting pairs of branch points of a 
hyperelliptic curve. The shape of the arcs is uniquely determined by a 
reality condition. For larger fi or for potentials with degenerate critical 
points, one expect the arcs to combine into graphs in the complex 
plane. It would be interesting to understand what kind of graphs can 
arise in this way. 
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